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Abstract. We investigate how the lensing potential can be measured tomographically with 
future galaxy surveys using their number counts. Such a measurement is an independent test 
of the standard ACDM framework and can be used to discern modified theories of gravity. 
We perform a Eisher matrix forecast based on galaxy angular-redshift power spectra, assum¬ 
ing specifications consistent with future photometric Euclid-like surveys and spectroscopic 
SKA-like surveys. Eor the Euclid-like survey we derive a htting formula for the magnifica¬ 
tion bias. Our analysis suggests that the cross correlation between different redshift bins is 
very sensitive to the lensing potential such that the survey can measure the amplitude of 
the lensing potential at the same level of precision as other standard ACDM cosmological 
parameters. 
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1 Introduction 

Cosmology has become a data driven science. After the amazing success story of the cosmic 
microwave background (CMB), see [1-5], we now also want to profit in an optimal way 
from present and future galaxy catalogs. Contrary to the CMB which comes from the two 
dimensional surface of last scattering, galaxy catalogs are three dimensional and therefore 
contain potentially more, richer information. 

In the large galaxy surveys planned at present one distinguishes spectroscopic surveys, 
which determine the redshift of galaxies very precisely and are used to determine the large 
scale structure (LSS), i.e., the clustering properties of galaxies, and photometric surveys 
which determine the redshift with less precision but which are optimized for galaxy shape 
measurements. Erom shape measurements one can then statistically infer the shear 7 and 
from it the lensing potential. In this paper we study to which extent the lensing potential 
can be obtained by simply looking at the angular correlation function of galaxies at different 
redshift without invoking shape measurements which are plagued by difficult systematics and 
intrinsic alignment [ 6 ]. 

The lensing potential is especially interesting since its relation to the matter density is 
determined by Einstein’s equations. The lensing potential ip is given by (see, e.g., [5]) 

'ip{n,z) = - df ^ ^ ^^ (<1> + T)(fn, To - f) (1.1) 

Jo r{z)r 
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where r(z) is the comoving distance of the source and tq is the present conformal time. 
We neglect spatial curvature, = 0, and we define the lensing potential as a function of 
redshift z and observer direction n. The variables <1> and 'h are the well known Bardeen 
potentials [5, 7]. By measuring both, the density fluctuation (5(n, z) and the lensing potential 
we can in principle test General Relativity on cosmological scales. 

When observing galaxies we measure their redshift and angular position. We can then 
determine the number of galaxies per solid angle and per redshift bin. The correlation 
function of these number counts within linear perturbation theory has been determined in 
Refs. [8, 9]. Apart from the galaxy number density and velocity it also depends on the 
convergence 

K = -^Ani/), ( 1 . 2 ) 

which affects the volume corresponding to a given observed solid angle and redshift bin. Here 
Aq is the angular Laplacian. 

Number counts have also been studied in [10-13] and relativistic expressions to second 
order have been derived in [14-18]. Their potential for cosmological parameter constraints has 
been analyzed in several papers [19-26]. In Ref [27] a code (CLASSgal) for fast computation 
is presented and forecasts for DES and Euclid-like catalogs are compared with the tradi¬ 
tional LSS analysis. In Ref. [28, 29] the capacity of SKA (the square kilometer array [30]) 
to determine primordial non-Gaussianity is analyzed using number counts. Eurthermore, 
relativistic contributions to the number counts can be isolated using special observational 
techniques [31-33]. In most of this paper we are concerned with the well known convergence 
term k which strictly speaking is also a “relativistic contribution”. However, due to the 
Laplacian, on small and intermediate scales, this term is much larger than the sub-leading 
velocity terms and the new “relativistic terms” derived in [8-10]. The amplitude of the latter 
is of the order of the gravitational potential, we therefore call them “potential terms”. We 
shall see that they are relevant only at very large angular scales and we shall neglect them for 
most of this work. Analyses of relativistic number counts in the context of modified gravity 
can be found in [34, 35]. 

Gosmic magnification /x = [(1 — k)^ — |7|^] ^ ~ 1 -|- 2ac [6 , 36, 37], has been detected 
by correlating background quasars and foreground galaxies, e.g., with the SDSS survey [38] 
or other catalogs [39-42]. Recent observations include, e.g., a detection of a redshift-depth 
enhancement of background galaxies magnified by foreground clusters using BOSS-Survey 
galaxies [43], and a measurement of the effects of lensing magnification on the detected 
number counts of background luminous red galaxies (LRGs) by foreground LRGs and clusters 
by [44] for the MegaZ (SDSS DR7) catalog. 

The analysis presented in this paper goes beyond this work. We constrain the lensing 
potential using the full tomographic redshift information of the lensing convergence, which 
allows us to probe the 3D information of a galaxy catalog. In certain situations the lensing 
term actually dominates the galaxy number count fluctuations and future surveys can be 
used to determine it. Our method is also very complementary to the usual determination of 
the lensing potential via shear measurements. It has different systematics and measures a 
somewhat different observable. 

This paper is structured as follows: in the next section we study number counts and 
determine the situations in which the lensing term dominates. In section 3 we present a Eisher 
matrix study for a photometric Euclid-like and a spectroscopic SKA-like survey. In section 4 
we conclude. Appendix A shows the survey specifications that we assume. In appendix B 
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we compute the magnification bias for Euclid. In appendix C we present some properties of 
spin weighted spherical harmonics used in the main text and appendix D outlines our Fisher 
matrix formalism. 


2 Determining the lensing potential with number counts 

In this section we recall the full expression determining galaxy number counts. We then 
analyze the different contributions to the power spectra and show how we can isolate the 
correlation of the density fluctuation with the lensing term. We also discuss the physical 
importance of a measurement of lensing convergence. 


2.1 Galaxy number counts 

Number counts of a given species of objects, e.g. galaxies is given by n( 2 :, n) = n{z)[l + 
A(z, n)] where n[z) is the mean galaxy density per redshift and per steradian at redshift z 
and [8, 9, 19, 28] 


A(n, z, miim) = b{z)D + — 




+ (2 - 5s) 


L ./0 ^ 


+ (/evo - ^)nv + 


/ 7Y 9 — 5s 


^ + drV+ / + 

Jo 


( 2 . 1 ) 


A dot denotes a derivative w.r.t. conformal time, Ti is the conformal Hubble parameter, V 
is the velocity potential for the peculiar velocity in longitudinal gauge, Vi = —diV, and D is 
the matter density fluctuation in comoving gauge while b{z) denotes the galaxy bias. More 
details are found in [19] and [27]. 

Denoting the limiting luminosity and magnitude of the survey by Liim and mum respec¬ 
tively, we have introduced the evolution bias, which captures the fact that new galaxies form 
and galaxies merge as the universe expands, hence their number density evolves not simply 
like (1 -|- z)^. 


din (a^n{z, L > Liira)) 

JevoyZL\iYn_) = 

a In a 


( 2 . 2 ) 


where n{z, L > L\\^) indicates the number density (per redshift and per steradian) of galaxies 
with luminosity above L\\^ and a is the cosmic scale factor. We also consider the magnifica¬ 
tion bias which takes into account that due to magnification less luminous galaxies still make 
it into our survey if they are in a region of high magnification and vice versa, 




dlogiQ n{z,L > Liim) 
dm 


miin 


(2.3) 


see appendix B for more details. 

The functions fevo{z, Tum) and s{z, miim) have to be determined from the catalog spec¬ 
ifications, either directly from observations or using simulations. As an example, in [28] the 
determinations of 5(z, Lum), /evo(2:, Tiim) = be{z,L\i„J) and 5s(z,miim) = 2(5(2, Lum) from 
simulations for the SKA survey of the 21cm line are discussed. Our modeling of these func¬ 
tions is described in appendices A and B. 
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The second term in the second square bracket of eq. (2.1) is the leasing term, 

A^{n,z) = —{2 — 5s)k . 


(2.4) 


Note that this term is specific to number counts. It comes from the change of the 
transverse surface area which has to be taken into account when relating number counts 
to density fluctuations. Due to the reciprocity relation this term cancels (at linear order 
in perturbations theory) in intensities like the 21cm intensity mapping [45] or the CMB 
anisotropies [5]. 

We can now expand A(n, z, miim) in spherical harmonics and compute its power spec¬ 
trum. Suppressing the limiting magnitude (or luminosity) mum we have 

^(n,^) = ^aemiz)Ye^(n) (2.5) 

i,m 

= 5(,ii5ynm'Ci{z,z'), ( 2 . 6 ) 

the Kronecker-deltas are, as usual, a consequence of statistical isotropy. 

In a true catalog we have to use redshift bins of finite thickness. We denote a normalized 
window function with width Az centered around the redshift Zi by Wi^z{z^ zi) and introduce 
the correlation power spectra 

C\^ = Ci{i,j) = I dzdz'WAziz, Zi)WAz{z', Zj)Ce{z, z') , (2.7) 

which has been thoroughly studied in literature [19, 46, 47]. In these observables, the dif¬ 
ferent terms of eq. (2.1) contribute differently. The first term, the standard density term 
A^(n, z) = bD, usually dominates. The standard redshift space distortion term in the 
Kaiser approximation 

A^^‘^{n,z) = n-^d^V , ( 2 . 8 ) 

is also very significant, especially for narrow window functions, Az < 0.01(1 -|- z). We stress 
that the usual Kaiser approximation for the power spectrum in Fourier space P{k) or the 
correlation function in configuration space also assumes that velocities of each galaxy pair are 
parallel, which implies n* = nj for each pair i,j, hence vanishing angular separation. Such a 
“flat sky approximation” is not necessary for angular correlations, thus including A’’®'^(n, z) 
in eq. (2.7) we also take into account so-called wide-angle effects [23, 24, 48-50]. We denote 
all the other velocity terms as “Doppler terms”, A'^(n, z) and all remaining terms apart from 
the leasing contribution as “potential terms”, A^(n, 2 ;), so that 

A = A^ + A"®^ + A^ + A^ + A^ with 

( TY 2 — "is \ 

w + 

2 — 5s 7” 

A^(n, z) = H-/ df (4> -7 T) -b (5s - 2)4> -b T -b 

+ (^ + ^ + + J'dfH + i)) 

The integrals in the first and second line of eq. (2.11) correspond to the Shapiro time-delay 
and integrated Sachs-Wolf effects, respectively. 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 
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Figure 1. Auto-correlations of the 4*^ bin (see figure 10 for an illustration of our binning). Due 
to the large width of the bins, redshift-space distortions are only important at large scales. On the 
other hand, lensing is the second important effect. Auto-correlations are dominated by the density 
term. Dashed lines show the sum of the auto-correlation of a given term plus its correlation with the 
density (we plot absolute values). 


We compute the power spectra in eq.(2.7) using the Class code^ [51]. We generalize 
version 2.3, already including the CLASSgal modifications [27], to take into account the 
redshift dependence of the bias parameters (appendix A). Furthermore, even if the expressions 
given above are valid at linear order in perturbation theory, we compute approximate non¬ 
linear spectra obtained using Halofit [52] as described in appendix D. 

2.2 The contributions to the number counts 

In the rest of this section, except when mentioned explicitly, we take as reference case the Eu¬ 
clid photometric specifications for 10 bins containing equal numbers of galaxies as described 
in appendix A. 

In figure 1 we show the diagonal spectra of the pure density term, Cf{i, i), the redshift 
space distortion term, C'^®^(f,i), the lensing term, the Doppler term Cf {i,i) and 

the potential term Clearly, the density term dominates, followed by the redshift 

space distortion and, on smaller scales, by the lensing term. Doppler and potential terms 
are significantly smaller. This figure is useful to compare the relative importance of different 
effects, but for the analysis of the full signal also cross correlations between the various terms 
must be considered. Since density is the dominant effect, in figure 1 we plot as dashed lines 
the spectrum of the correlation of a given term X plus its cross-correlation with the density, 
schematically \ {X{zi)X{zj)) + {X{zi)5{zj)) + {5{zi)X{zj))\. As we shall see below, the lensing 
term is always dominated by the density-lensing cross correlation rather than the lensing¬ 
lensing correlation itself. We have verified that in the case of the spectroscopic SKA surveys, 
thanks to the good redshift determination, galaxies are not spread over Gaussian bins and the 
relative importance of local terms (“D”, “rsd”, “V”) increases compared to integrated effects 
(“L”,“P”). See also [19] for a comparison between photometric and spectroscopic surveys. 

^http://class-code.net/ 
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Bins 1-10 




Bins 1-4 Bins 1-4, SKA 




Figure 2. Top left panel: the auto-correlation of the 1*^ bin (see figure 10 for the binning), top right 
panel: the cross correlation of the bins 1-10 and bottom left panel: the cross correlations of bins 1-4 
are shown for an Euclid-like survey. For comparison, in the bottom right panel the correlation of bins 
1-4 are shown also for SKA. The solid (black) line “full” contains the correlation of density, redshift 
space distortion (rsd) and lensing; the dashed (blue) line “only local” includes density and rsd only; 
the short-dashed (violet) line “local-lensing” is the correlation of the local terms in bin i with the 
lensing term in bin j, with mean redshift Zi < Zj. 


In our Fisher matrix analysis in Section 3 we neglect the potential terms since the 
likelihood is strongly dominated by £ > 20 where they do not contribute appreciably. 

2.3 Measuring the cross correlation {D{zi)K{zj)) 

The situation changes drastically if we consider the cross correlation spectra. In figure 2 we 
show them for the density, redshift space distortion and lensing terms neglecting the sub¬ 
leading Doppler and potential terms. The “full” case contains the correlation of all these 
terms, compared to the case where we consider only “local” terms (density and redshift-space 
distortions) and the correlation of local terms in bin i with lensing term in bin j with mean 
redshift Zi < zj (“local-lensing” case). Note that when fixing Zi and increasing Zj, while the 
density and redshift space distortion terms decay, the lensing term remains substantial. We 
have also checked that the Doppler and potential terms remain small for £> 20, and since 
larger scales are dominated by cosmic variance we neglect these effects to perform forecasts, 
letting a detailed study of their significance as a separate work. The figure shows that bin 
auto-correlations are dominated by local terms. However, lensing convergence dominates the 
cross spectra. The correlation 1-10 is entirely given by the density-lensing correlation (short- 
dashed). For the 1-4 correlation, the density auto-correlation is still important because 
of the significant overlap of the wide bins 1 and 4 due to the poor photometric redshift 


- 6 - 


























determination. For comparison we also show the correlation of bins 1-4 for SKA: given 
the spectroscopic redshift precision, this correlation is already dominated by the leasing 
contribution. 

The largest contribution comes in particular from correlating the density transfer func¬ 
tions at a lower redshift bin i with the integrated leasing transfer function at higher redshift 
bin j, with i < j. This corresponds to a measurement of leasing of background sources 
by foreground sources. The inverse case i > j gives a negligible contribution. This shows, 
that for sufficiently distant bins or bins with negligible overlap like the 1-4 bin of SKA, the 
galaxy number counts actually measure the following quantity averaged over the respective 
bins with the window functions WA.z{z,Zi) and WAz{z',Zj) 

-b{z){2-5s{z')){D{z,n)K{z',ii')) = + l)Ci{z, z')Pe{n ■ n). (2.12) 

i 

Here Pi is the Legendre polynomial of degree i. All other terms contribute very little to the 
number count fluctuations if r{z') — r{z) > Vc- Here Vc is a typical galaxy correlation scale of 
order 150/i“^ Mpc, and for a mean redshift ^ ~ 1, r{z') — r(z) > (1.76 x Hq)~^{z' — z), this 
corresponds to z' — z > 0.09. For well determined redshifts, like e.g. SKA this is satisfied 
for all not neighboring bins in our binning. In the case of photometric redshifts the overlap 
of the different bins is considerable and the local terms still contribute significantly also to 
bins with well separated mean redshift as can be seen in the lower left panel of Fig. 2. We 
have verified that for the cross-correlations shown in figure 2, rsd’s contribute only to low 
multipoles £ < 10 dominated by cosmic variance. However, for narrower bins, rsd can be 
relevant and should be added to the density term in eq. (2.12). 

For a given (reconstructed) bias b{z) and magnihcation bias s(z), the cross correlation 
spectra of number counts allow us to determine the power spectrum of {D{zi, n)K{zj,n')) for 
sufficiently well separated redshifts, Zj — Zi > 0.1. In the next section we shall show in a 
simple example how bin cross-correlations can be used to constrain modihcations of General 
Relativity on cosmological scales. 

Note also the sign change between the top and the lower right panels in Fig 2. This 
is due to the sign change in 2 — 5s(z') which happens at s{z') = 0.4. The cross correlation 
spectrum of {D{z, n)K{z', n')) is always positive so that Ci{z, z') given in eq. (2.12) is negative 
for s{z') < 0.4 and positive for s{z') > 0.4. In our two examples, 2 — 5s(2;) becomes negative 
for SKA in the 4th bin, while its mean remains positive for Euclid until the 8th bin. The 
dependence of s{z') on the limiting magnitude can actually be employed to monitor this term 
by varying m\\^, see appendix B. 

Let us hnally stress that the oscillations visible in the spectra are not due to numerical 
precision but represent the Baryon Acoustic Oscillations visible mainly in the cross correlation 
spectra of the local terms. 

In hgure 3 we show separately the contributions (2 — 5si)(2 — 5sj)(KiKj} and —bi{2 — 
5sj){DiKj) to the cross-correlation of the redshift bins 1 and 10. The plot shows that the main 
contribution is the correlation of the density field at lower redshift with the lensing at higher 
redshifts, i.e., Zi < Zj. This is nearly 2 orders of magnitude larger than the lensing-lensing 
correlation, and about 6 orders of magnitude larger than the density-lensing correlation for 
Zi > Zj. In this last term numerical oscillations are clearly visible, but since the amplitude is 
negligible, they are not relevant for our Fisher matrix analysis in the next section. The results 
in figure 3 are expected, since lensing at Zj > Zi is caused by all galaxies up to redshift Zj 
hence also by the galaxies at redshift Zi which are represented by Di. This explains also why 
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Bins 1-10 



1 10 100 1000 


- -b-i(2-5sio)<Di/Cio> 

- _(2-5si)(2-5sio)<KiKio> 

. -(2-5si)bio</CiDio> 


/ 


Figure 3. The contribution from the density and convergence terms has been isolated for Euclid bins 
1-10 cross-correlation. Note that since (2 — 5si)(2 — 5sj){KiKj) is negative for the chosen values of 
magnification bias, we plotted it with opposite sign. Due to the sign change of (2 —5s(z)) at z « 1, the 
lensing-lensing correlation is negative, while the density-lensing correlation is positive. The numerical 
oscillations in (kiDiq) are not relevant for the analysis, as the amplitude of this term is negligible. 


(DiKio) ^ (kiDiq). Since the density term is by far the most relevant one, we also expect 
(DiKio) ^ (ki^io). Again, the sign of the correlation, negative for (2 — 5si)(2 — 5sj){KiKj) 
and positive for —bi{2 — 5sj){DiKj), is also determined by the sign of the factor (2 — 65 ( 2 ;)) 
which changes at z ~ 1. We have also verihed that our results for the convergence spectrum 
are consistent with, e.g., [53]. 

In figure 4 we also include the noise amplitudes (shaded regions) for the photometric 
Euclid survey with 5 and 10 Gaussian bins containing equal numbers of galaxies, see ap¬ 
pendix A. In this plot, for illustrative purposes, we estimate the error of the signal Cf as its 
variance by setting (ij) = (pq) in eq. (D.l): 

(2.13) 

where /sky is the covered fraction of sky and are the observable spectra including 

errors, see Appendinx D for more details. Errors are dominated by cosmic variance on large 
scales, and by theoretical errors Ei and shot-noise (see eq. (D.2)) on small scales. 

For the first bin auto-correlation, both conhgurations show good signal-to-noise. In 
the case of 5 bins, the signal is slightly lower. In this case the bins are wider, hence the 
signal is spread out over a larger redshift range and local terms (not integrated along the line 
of sight) dominating auto-correlations are somewhat less important. The cross-correlation 
between the first and last bins show that, despite the large redshift separation, the signal is 
still substantial both for 5 and 10 bins because of photometric overlap but also thanks to the 
contribution of lensing convergence as discussed above. The 5 bin case has the advantage 
that each bin has a lower shot-noise, which dominates the error at higher multipoles. 

In figure 5 we show the redshift integrated spectra, i.e., we choose only one bin covering 
the full range 0 < 2 ; < 2 of the survey. Given that the single redshift bin is now much 
larger than the photometric errors, we use a tophat window function. Furthermore, we now 
consider the full redshift dependence of galaxy and magnification bias within the window 


^ij 


/sky(2£+l) L 


^ D ^ C 


+ c: 


-iObs,ij 
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Figure 4. The total signal (solid line) together with the noise (gray band) are shown. The noise is 
dominated by cosmic variance at low multipoles and by shot-noise and theoretical errors at high 
multipoles (see appendix D for details). Upper and bottom panels consider a division into 10 and 5 
equal shot-noise redshift bins, respectively. Bin i-j correlations are indicated inside the plot. Thanks 
to lensing convergence, also cross-correlations have an important signal-to-noise for (. > 100 — 200. 


0<z<2 



D-D 

rsd-rsd 

L-L 

P-P 

V-V 


Figure 5. The redshift integrated spectra. Compared to the 10 bin configuration, local terms are 
significantly reduced, while integrated terms are nearly unchanged so they become more relevant. The 
notation follows figure 1. 
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function, unlike the cases for 5 and 10 bins where these parameters are set to their values at 
the mean redshift of each bin (see appendix A). Comparing the amplitude with the one in 
figure 1 shows that local terms are significantly reduced when integrating over a wide window 
function. On the other hand, terms integrated along the line of side change only little, so 
that their relative importance increases. Only the density term is still larger than the lensing 
term, while rsd and Doppler are subleading. The potential terms now even dominate the 
signal at ^ = 2 and i = 3 and are the second most important term for i < 10, whereas the 
lensing term is second for smaller scales. As for figure 1 , we stress that the full signal is not 
given only by the auto-correlation of each term (solid lines), but we have to add all cross¬ 
correlations between different terms. Dashed lines include also cross-correlations of each term 
with the dominant density term. Neglecting lensing would underestimated the total signal 
by 5-10%. On large scales, the potential terms become very important and can mimic a 
significant /NL-contribution to the bias, see [28, 54] for a detailed study of this effect. 


3 Lensing constraints for modified gravity from tomographic number counts 


In this section we show that the number counts can be used to determine the lensing spectrum. 
For this we modify the lensing spectrum by one simple parameter, 

V^(zi, n) = /3 n) . (3.1) 

In the standard ACDM model we have (3 = 1. In the literature [55, 56] one often finds the 
so called “slip parameter” -q and the “clustering parameter” Q given by 

$ = = dvrGa^gD. (3.2) 


Here and T are the Bardeen potentials and both r] and Q can in principle depend on time 
and on wave number. Inserting this in eq. (1.1) we obtain that with these modifications, for 
constant values of Q and ij and neglecting the anisotropic stress from neutrinos 

'*/’(^i, n) = ^Q(l -h n), (3.3) 

so that (3 = + r]~^). A deviation of this value from unity requires a modification of 

the standard ACDM cosmology. This can in principle be achieved, e.g., with clustering dark 
energy, but in most cases requires a modification of General Relativity [55]. Generally, also 
the parameter (3 may show a redshift and scale dependence that we neglect here for simplicity. 

We recall that the lensing convergence k is a complementary probe to the lensing shear. 
Convergence measures the trace of the deformation matrix [57] 


^ ^ d{0s,V’s) ^ / 1 - K - 7i -72 \ 

d{0o,^o) V -72 + 

while shear is defined as the traceless part. Setting 


7 = 


fli 72 ^ 
V 72 -7i J 


(3.4) 


(3.5) 


The tensor field 7 has helicity +2 with amplitude 71 -|- ^72 which is also denoted by 7 . 
Assuming a purely scalar gravitational field with lensing potential i(> given in eq. (1.1), so 
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that Aij = 5ij + we have 


(3.6) 

(3.7) 


K = = + 0* 0)^ 

7 = 7i + ^72 = 0^'0 ■ 


Here 0 is the spin raising and 0* the spin lowering operator, see appendix C and [5] 
for details. Expanding k and ijj in scalar spherical harmonics and 7 in spin-2 spherical 
harmonics one finds the following relations for their power spectra (the derivation is given in 
appendix C): 


C7 = 






- 

I'-'i — 


4 ' 


(3.8) 


Weak leasing experiments measure shear and hence 2 CJ while the proposed leasing tomog¬ 
raphy measures C^. The relations (3.8) therefore provide a welcome consistency check. 

We determine error bars on the leasing parameter 0 from a Fisher matrix analysis. We 
study the dependence of our results on the number of bins for both, a Euclid-like and an 
SKA-like survey. We choose the Planck 2015 results [2] as our fiducial parameters around 
which we compute the Fisher matrix. These are the following: cub = = 0.02225, Wcdm = 

= 0.1198, Us = 0.9645, ln(10^°As) = 3.094, Hq = 67.27km/s/Mpc, rrii, = 0.06 eV 
and 0 = 1, corresponding to the baryon and CDM density parameters, the primordial scalar 
spectral index, the amplitude of primordial curvature perturbations, the Hubble parameter, 
the mass of one single neutrino species (neglecting the other neutrino masses) and our leasing 
parameter 0. In joint constraints of a subset of parameters, we marginalize over the others. 

We divide the galaxies into A^bin = 1, 5,10 redshift bins with the specifications described 
in appendix A. We assume constant galaxy and magnification bias in each bin, the values 
being determined by the mean redshifts. Only for A^bin = 1 these parameters are fully evolved 
within the bin. 

We neglect the contribution of potential terms to the number counts, eq. (2.11), relying 
on the fact that, as shown in figures 1 and 5, they are only relevant at very large scales £ < 20. 
These terms include effects integrated along the line of sight as leasing, and neglecting them 
may bring systematic errors on the determination of the leasing parameter 0. On the other 
hand, effects integrated along the line of sight are computationally costly and considering 
them, e.g., in a Markov chain Monte Carlo may be challenging. For our purposes, actual 
observations can simply neglect scales ^ < 20 which, due to cosmic variance, contribute 
anyway very little to the constraining power. 


3.1 Euclid forecasts 

Figure 6 shows the forecasted 1-cr contours for a Euclid-like survey. Comparing the cases 
A^bin = 5 (red contours) and A^bin = 10 (black contours), we conclude that a larger number 
of bins clearly improves the constraints as we expect since we add more redshift information. 
Furthermore, considering only bin auto-correlations and neglecting bin cross-correlations 
(dashed contours) gives of course worse constraints than when including all correlations (solid 
contours). This is particularly relevant for the lensing parameter 0 since, as shown in figure 2, 
auto-correlations are only weakly sensitive to lensing, so the constraining power comes mainly 
from cross-bin correlations where the density-lensing correlation is the dominant contribution. 

In figure 7 we show the 1-dimensional marginalized constraints. In this case we also 
consider Abin = 1- Given the large width of the bin relative to photometric errors, in this case 
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Figure 6. Euclid 1-a (68.3%) 2d likelihood contours and the Id likelihood functions for 5 (red) and 
10 (black) redshift bins, including all bin auto- and cross-correlations (solid) or only auto-correlations 
(solid). For the ln(10^°As) - contour, the black full correlations (solid) 10-bin contour is slightly 
larger than the corresponding 5-bin contour. We have checked that this is a precision issue which is 
not relevant for the conclusions of this work. Note, however, that the leasing parameter /3 is mostly 
constrained by cross-correlations. It is the only parameter for which the 10-bin auto-correlation result 
(black dashed) is worse than the 5-bin cross correlation (red, solid). Also, increasing Win from 5 to 
10 only slightly improves the precision of /3 while it significantly improves the constraints on the other 
parameters. 


we assume a tophat window function covering 0 < z < 2. The case Win = 1 is characterized 
by errors about 1-2 orders of magnitude larger than for Win = 10. For the lensing parameter 
/3 the error increases even by 3 orders of magnitude, confirming that most of its constraining 
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Figure 7. Marginalized errors (68.3%) for Euclid, normalized to the fiducial value of a given 
parameter Aq,. On the cc-axes the numbers corresponds to iVbin and the labels indicate the situation 
in which all bin correlations are taken into account versus considering only auto correlations. The red 
line ln(10^°As) is mainly sensitive to Agb^, where &o is the bias at z = 0. Since we have set Bq = 1, 
marginalization over bo would substantially enhance the error of Ag and, consequently, on Ug. Note 
that the precision of /? is comparable to the one of the other parameters. 


power comes from cross-correlations of different bins. For A^bin = 5,10, the forecasts suggest 
0(1%) constraints on /3. However, because of the limitations of a Fisher analysis, rather than 
relying on the absolute constraint itself, we compare it to those of other standard ACDM 
parameters to conclude that we expect competitive measurement of /3. We have verified that 
this still holds for the more pessimistic assumption on nonlinear scales setting ^max = 1000 
instead of 2000 (see appendix D), for which 1-dimensional error bars of each parameter 
typically increase by a factor 2-3. 

As explained in appendix A, photometric redshift uncertainties can be modeled by Gaus¬ 
sian bins (centered at the estimated true redshifts) with variance determined by the typical 
photometric error. This reduces the signal, which is integrated over a broader redshift range. 
In general, Gaussian redshift bins overlap and, as shown in figure 2, the relative contribution 
from the lensing signal can be reduced for some cross-correlations, compared to the case 
of non-overlapping tophat bins (e.g., the spectroscopic SKA). On the other hand, lensing 
increases coherently with the width of the redshift bin while the density signal decreases. 
Hence, distant correlations with small overlap are still important enough to lead to a clear 
improvement in the determination of [3 when including all cross-correlations, as shown in 
figure 7. In the next section we provide a similar study for the case of the spectroscopic SKA 
with tophat bins. We shall see that in both cases, including cross-correlations reduces the 
error in /? by a factor of about 3 for the case of 5 bins and by a factor of about 2 for 10 bins. 
There is no significant difference in this for Gaussian or tophat bins. 

It is interesting to compare the present forecasts to those from standard galaxy clus¬ 
tering and weak lensing analysis. We refer to section 1.8 of [58] and to [59]. Euclid is 
expected to measure the main standard cosmological parameters to percent or sub-percent 
level. Our constraints on the lensing parameter /3 can be compared to those on modified 
gravity parameters for the simplest models, as outlined in eq. (3.3). Modified gravity con¬ 
straints are expected from analysis of redshift-space distortion through the measurement of 
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Figure 8. SKA 1-a (68.3%) 2d likelihood contours and the Id likelihood functions for 5 (red) and 
10 (black) redshift bins, including all bin auto- and cross-correlations (solid) or only auto-correlations 
(dashed). Also here some of the black full correlations (solid) 10-bin contours are slightly larger than 
the corresponding 5-bin contours due to precision. Again, for the lensing parameter /3 which is mostly 
constrained by cross-correlations, the 10-bin auto-correlation result (black dashed) is worse than the 
5-bin cross correlation (red, solid). 


power spectra or correlation functions in Fourier space. These observations perform better 
when spectroscopic redshifts are provided and constrain the matter growth factor index 70 - 
defined by 73 ^ = ^ where G{a) is the linear growth function. The relation of 

7 g to the parameters in eq. (3.2) can be found in, e.g., [60]. The index 7 ^ is constrained 
by few percent for the simplest models, but the constraints weaken when marginalizing over 
galaxy bias, including a redshift dependence yc = 7 g(^) or when marginalizing over several 
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Figure 9. Marginalized constraints (68.3%) for SKA, normalized to the fiducial value of a given 
parameter Aq,. On the x-axes the numbers corresponds to A^bin and the labels indicate the situation in 
which all bin correlations (“all”) or only bin auto-correlations (“auto”) are included. As for figure (7), 
marginalization over Bq would substantially increase errors especially on Ag and, consequently, on ng. 


modified gravity parameters. The constraints can easily reach 0(10%). We expect a similar 
weakening of the constraints for /3 when introducing, e.g., a possible redshift dependence of 
13. Nevertheless, also in the simplest case, when /3 = 1 is inconsistent with data, this already 
implies a deviation from the standard theory of gravity. Therefore, this constraint can be 
viewed as a consistency check of GR. 

Weak leasing constraints via two-point tomographic cosmic shear measurements are 
also very interesting for the present analysis. The Euclid shear power spectrum is expected 
to be recovered to sub-percent accuracy over signal-dominated scales, giving an integrated 
accuracy (over all scales i < iraax ~ 5000) of 10“'^ [59]. The main difficulties of shear 
measurements are the modeling of non-linear scales and of intrinsic alignment. On the other 
hand, our measurements of the magnification k depend significantly on the modeling of galaxy 
and magnification bias, but do not suffer from intrinsic alignment. Therefore, checking the 
relation (3.8) is also an excellent test for systematics. Finally, we note that while shear 7 is 
an observable, convergence k is not gauge invariant [ 8 , 61], but figures 1 and 5 prove that 
the terms involved in gauge transformations are subleading in our configurations. 

3.2 SKA forecasts 

Figure 8 shows the forecasted 1-a contours for the SKA survey. Both Abm = 5 (red contours) 
and Abin = 10 (black contours) show that including all bin correlations (solid lines) improves 
the constraints obtained when considering only bin auto-correlations (dashed lines). Fur¬ 
thermore, in this case with very precise redshift determination, redshift bins do not overlap 
and are more weakly correlated compared to the photometric Euclid case. In this case, the 
residual cross-correlation is nearly entirely due to integrated effects. 

Comparing the Abin = 5,10 contours (red lines to black lines) shows that increasing 
the number of bins always improves or gives similar constraints. By running parts of our 
simulations at higher precision, we have checked that the slight exceptions to this in some 
of the 2d contours are due to precision issues. The improvement with the number of bins 
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is also shown in figure 9, where we plot 1-dimensional errors. Furthermore, hgures 8 and 9 
show that the j3 parameter improves significantly when including all bin cross-correlations. 

Similarly to a photometric Euclid-like survey, also the analysis for SKA suggests a 
constraining power on the lensing parameter similar to the errors on the other standard 
ACDM parameters. We have also verified that, as for Euclid, this relative comparison still 
holds for the more pessimistic assumption on the non-linearity scales setting ^max = 1000 
instead of 2000 (see appendix D), the error on each parameter again increases by a factor 
~ 2 in this case. 

Given that for A^bin = 1 (0.1 < z < 2) the differences between Euclid photometric and 
SKA spectroscopic redshifts are no longer important, we expect similar constraints from SKA 
as those obtained in figure 7 for Euclid. 

4 Conclusions 

We have investigated how the lensing potential can be constrained tomographically using the 
galaxy number counts statistics in future surveys. We use the dependence of the observable 
number counts on the lensing convergence k. Redshift bin auto-correlations are dominated 
by local terms like density and redshift-space distortions, which give small contributions to 
redshift bin cross-correlations. Redshift bin cross-correlations separated by large radial dis¬ 
tances mainly come from contributions integrated along the line of sight, which are dominated 
by the lensing convergence. The number count cross correlation spectra with little overlap 
provide an excellent measure of the power spectrum of {D{zi,n)K{zj,n')) for Zj > zi with 
Ar = r{zj) — r{zi) > 150/i“^Mpc, hence Az = Zj — Zi> 150 Mpc/i“^F7(z) « 0.09 for z ~ 1. 

Eurthermore, for very deep redshift bins (0 < z < 2), lensing is a significant contribu¬ 
tion also for the bin auto-correlation, confirming what has already been found in previous 
work [19]. However, parameter constraints are in this case nearly 2 orders of magnitude worse 
compared to the subdivision into Abin = 5,10 which makes use of the redshift information. 

The cross-correlations term dominating number counts for well separated redshift bins, 
given by the density-lensing correlation {D(zi)K{zj)), corresponds to the lensing of back¬ 
ground galaxies by foreground galaxies, and it is expected consistently with observations 
(e.g. [38, 43]). Lensing magnification studies usually work with two different populations, 
one providing the foreground and the other the background galaxies. Here we go beyond this 
approach and consider a fully tomographic extraction of the lensing signal as a function of 
redshift. Since the contribution {K{zi)D{zj)) with Zi < Zj is negligible, there is no danger of 
contamination. Given D{zi) from the auto-correlation spectrum, the cross-correlations term 
provides an excellent measure of the lensing potential via 

{D{Zi)K{Zj)) 

' 

As an application in this paper we have introduced the parameter (3 through eq. (3.1). 
The standard AGDM model expects (3 = 1 and a deviation from this value may indicate 
a modification of General Relativity. In any case, constraints on f3 provide a consistency 
test of the cosmological standard AGDM model, which can be considered together with 
the measurements of other parameters like, e.g., the growth factor or the equation of state 
[58] which would indicate deviations from AGDM. The proposed method provides also an 
independent measurement of the lensing potential, alternative to shear observations [6] and 
it is therefore a very useful consistency check. 
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We have performed a Fisher matrix forecast based on a tomographic analysis with 
angular-redshift power spectra. Nonlinear scales are included with a theoretical error which 
we model as in [62], see Appendix D. We expect that both photometric Euclid-like and 
spectroscopic SKA-like surveys will produce constraints on the lensing parameter j3 that are 
competitive with those on the other standard ACDM parameters, provided that all redshift 
bin cross-correlations are taken into account. This result has been obtained by considering 
realistic specifications, in particular we have derived a redshift-dependent fitting formula for 
the magnification bias, eq. (A.5), consistent with photometric Euclid-like surveys. 

As a simple continuation of our analysis one can fit the value of /3 in each redshift 
bin. This corresponds to a measurement of the lensing potential tpizj). However, given the 
limitations of a Eisher forecast, a definitive estimation of errors on the lensing parameter /? 
or I3{zi) should rely on a Markov chain Monte Carlo forecast. In such a forecast one would 
want to consider a specific experiment with its predicted likelihoods in more detail and also 
marginalize over the necessary nuisance parameters. To aim at percent level constraints, our 
approximate treatment of non-linear scales (i.e., the rescaling of linear transfer functions by 
the Halofit algorithm) should be also improved. Integrated effects other than lensing have 
been neglected in the present Eisher analysis, but to avoid systematic errors on the lensing 
parameter f3 from actual observations, they should either be consistently included, or the 
large scales £ < 20 where they are significant must be excluded from the analysis (note, 
however, that on this scales cosmic variance is substantial so that they do not contribute 
significantly to the final result). 

Eurthermore, it has been found that size bias, i.e., the fact that only galaxies large 
enough to be detected as extended sources are included in catalogs, may be even more 
relevant than magnification bias [63] and should be included to avoid systematic effects. 

We also consider a promising future possibility to change s(^) by varying the limiting 
magnitude of the survey and in this way to enhance or reduce the lensing signal. 
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Figure 10. Euclid (left) and SKA (right) galaxy density distribution (black line) with a division into 
10 bins containing the same number of galaxies. 


A Surveys specifications 
A.l Euclid 

Following [58, 59], we consider Euclid photometric specifications and approximate the number 
of galaxies per redshift and per steradian, the galaxy density, the covered sky fraction, the 
galaxy bias and magnification bias as 


ai\ o o 

—— = 3.5 X lO^z^exp 

dzdil 

for 0 < z < 2.0 , 
d = 30 arcmin”^ , 

/sky = 0.375 , 
b{z) = boVTTz , 
s{z) = So + SlZ + S 2 z‘^ + s^z^ . 



(A.l) 


(A.2) 

(A.3) 

(A.4) 

(A.5) 


where zo = Zmean/1-412 and the median redshift is Zmean = 0.9. We set 6 o = 1 throughout this 
work, but it is straight forward to include it or to marginalize over it. The magnification bias 
is computed in appendix B and the coefficients are sq = 0.1194, si = 0.2122, S 2 = —0.0671 
and S 3 = 0.1031. The full redshift and limiting magnitude dependence is given in eq. (B.9). 
Photometric redshift errors 6z = 0.05(1 + 2 ;) allow us to choose up to Abin = 10 bins such that 
the zth-bin width is A.Zi > 26z- Figure 10 shows the division into 10 bins containing the same 
number of galaxies, and where redshift uncertainties are taken into account by modeling the 
bins as Gaussian with standard deviation Azij2. For simplicity, also for Abin = 5 we assume 
Gaussian bins with standard deviation Azij2, even if in this case the bins are larger than 
photometric resolution (but still comparable). For numerical convenience we set the lower 
redshift bound to 2 ; = 0.1; this affects our results by a negligible amount. Figure 11 shows 
the redshift dependence of galaxy and magnification bias. Except when specified differently, 
we assume constant galaxy and magnification bias in each bin, the values being determined 
by the mean redshifts. 
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Figure 11. Galaxy bias b{z) and magnification bias s{z) for Euclid and SKA. Euclid’s magnification 
bias is computed for mu^ = 24.5. 


For simplicity, to estimate the evolution bias fem (that in principle can be computed 
from the luminosity function), we assume to observe all the galaxies in the windows and 
insert dN/dz/dVl in eq. (2.2) to compute it. Given the high success rate of the survey (larger 
than 50%) and the fact that /evo only appears in the subleading Doppler and potential terms, 
this approximation does not affect our results. 


A.2 SKA 

For SKA2 we use the specifications from [28, 64] 


^=(1^) 10«.«exp(-«^) (A.6) 

for 0.1 < z < 2.0 , 

/sky = 0.73 , (A.7) 

6 (z) = C 4 exp (csz) , (A.8) 

s{z) = ce + C 7 exp{-csz) , (A.9) 

where ci = 6.7767, C 2 = 2.1757, cs = 6.6874, C 4 = 0.5887, cs = 0.8130, ce = 0.9329, 
cj = —1.5621, Cs = 2.4377. dN/dz/dO, is the number of galaxies per redshift and per 
steradian. The magnification bias is computed from [28] considering that the specifications 
described above are consistent with the 5/rJy sensitivity. Given the spectroscopic redshift 
determination, we use tophat redshift bins. In principle order ~ 10^ bins can be considered. 
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but for comparison with the Euclid photometric case we limit ourselves to 10 bins, see 
figure 10. For an analysis of constraints improvement with the number of bins up to the limit 
set by shot-noise, see [65]. Also in this case, except when specified differently, we assume 
constant galaxy and magnification bias in each bin. 

B Euclid magnification bias 

In this appendix we derive a fitting formula for the magnification bias of the photometric 
Euclid survey, relying on estimates of the luminosity function. The fit is useful to perform 
forecasts for Euclid-like surveys and it is not meant to be employed for precise computations. 
For more involved and more accurate estimation in the optical bands see, e.g., [53]. 

The Schechter luminosity function is parametrized in terms of absolute magnitudes M 
as 


(l){M,z)dM = 0.41n(10)(/>*(z) 


100'4(^*(-)-^)r^'exp 




dM 


(B.l) 


where M*{z) is a characteristic magnitude, (t)*{z) is a number density and a is a constant 
faint-end slope. The Euclid photometric survey will observe in the visible and near-infrared 
bands in the range 550 — 900nm [59]. We use the results of [66] to model the redshift 
dependence of the luminosity function. We take the i' filter at ~ 770nm as reference. The 
Schechter parameters are given by (see Table 9, Case 3 of [66]): 


M*{z) 

— Mq -|- Hi ln(l -|- z) , 

(B.2) 

<P*iz) 

= 4>oi^ + ^)“^ ) 

(B.3) 

a 

= -1.33 , 

(B.4) 


where oi = —0.85, 02 = —0.66, Mg = —21.97 and = 0.0034 Mpc“^. 

The cumulative number density of galaxies n{z) = dN{z, M < M\\^)/dz/dO- with mag¬ 
nitude lower than M\i^ is: 


n(M < Miim) = 


'■Mim 


(j){M,z)dM . 


The magnification bias is defined as in eq. (2.3): 

Slogigh 


S(2,miim) = 


dm 


miin 


where m is the apparent magnitude related to the absolute M by 

'diiz) 


M = m — 5 log;^g 


10 Pc 


-Kiz), 


(B.5) 


(B.6) 


(B.7) 


where di is the luminosity distance. We also include iC-corrections due to the fact that we do 
not observe bolometric magnitudes (integrated over all frequencies), hence a fixed observing 
band dims with redshift. For a (x spectrum we have 


K{z) = 2.5(7 - 1 ) logio(l -b z) . 


(B.8) 










Figure 12. Euclid magnification bias at the mean redshift ^mean = 0.9, as a function of the limiting 
magnitude. 


Typical values in the B (ultraviolet) and K (infrared) bands are 7 ~ 4 and 7 ~ —1.5, 
respectively [67]. Comparing to the itT-corrections estimated in, e.g., [68], we deduce 7 « 2.7 
for our i! band. 

Given that, at fixed redshift, derivatives w.r.t. m coincide with derivatives w.r.t. M, 
we can write 


s{z,miiyn) 


d logio ^ _ 1 (/>(Miim, 2 :) 

dM In(lO) h(M < Mu^) ’ 


(B.9) 


Note however, that for fixed mum, Tfum depends not only on 2; but also on cosmological 
parameters via diiz). Furthermore, we have taken into account that di also contains fluc¬ 
tuations [27, 69]. 

Figure 12 shows the dependence of the magnification bias on the limiting magnitude 
miim at the mean redshift Zmean = 0.9. As expected, higher mum correspond to lower values 
of magnification bias. However, for large values of m\\^ magnification bias goes towards 
a constant value. If the faint-end of the luminosity function would have been flat (slope 
a = —1), then s —)■ 0 as m\\^ —)■ 00. 

In figure 11 we also show the magnification bias for the Euclid limiting magnitude 

= 24.5 at a function of redshift. This is well approximated by the fit given in eq. (A. 5). 
The values s{z = 0.5) = 0.22 and s{z = 1.0) = 0.37 are consistent with the results obtained 
in [70]. We stress that the critical value s = 2/5 for which leasing convergence vanishes, is not 
an issue for our analysis thanks to the redshift dependence of ^(z). One could actually use 
different cuts in magnitude that would lead to an optimized s{z) for our analysis. Especially, 
a slight increase of the limiting magnitude to e.g. m\\^ = 25 would substantially enhance the 
leasing contribution. 
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C Spin weighted spherical harmonic decomposition of the lens map 

An arbitrary tensor field on the sphere can be decomposed into its spin components. A 
pure spin s tensor held is a symmetric completely traceless rank s tensor and always has 
two helicity states, namely helicity ±s. These can be expanded in terms of spin weighted 
spherical harmonics The normal spherical harmonics Vim are the ^Yim- See [5] for 

more details. We therefore can expand the scalar functions n and and the helicity +2 held 
7 as 


i) 

W 

II 

(C.l) 


ira 


K 

— 0.£m^lrn 

(C.2) 


im 


7 


(C.3) 


Im 


Applying the spin raising operator twice on V’ one obtains the helicity +2 component of the 
Jacobian Aij of the lens map, 

7 = -2^V- (C-4) 

In terms of spin raising and spin lowering operators, the Laplacian is given by [5] 

= + $* =-2k. (C.5) 


Expressed in terms of the polar basis (??, (p) on the sphere the spin raising and spin 
lowering operators are given [5] $ = —\/2V_ and $* = —\/2V+, where V± denotes the 
covariant derivative in the directions of the helicity vectors e± given by 


y/2 Sind 


We now use the following properties of spin weighted spherical harmonics [5]: 


0 is^em) — \/{(■ — s){(. + s + 1) 

0* (gYem) = -V{e + s){£-s + l) g_,Y,^. 

This also implies 

AY^rn = \{ 00 * Y 0 * $)Ytm = -£{£ + l)Yim 

With this we hnd 




7 _ 

2^£m o 


(^ + 2 )! 

7 - 2 )! 


a 


£m ■ 


(C.6) 

(C.7) 


(C.8) 


Using the dehnition of the expectation values, Cf = (|a^P), one obtains eq. (3.8). 
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D The Fisher matrix 


We introduce the covariance matrix as in 


[46] 


^ohs^ip^ohs,jq , ^ohs^iq^ohs^jp 


AkyA£(2^+l) 


(D.l) 


where A£ is the bin width in multipole space. We choose IS.!!. = 2//sky (rounded to the closest 
integer value), so that the covariance matrix is approximately block-diagonal, as empirically 
demonstrated in [71]. The theoretical error on non-linear scales and shot-noise are taken into 
account as 

= + + (D.2) 

where is a theoretical error on power spectra, and M is the number of galaxy per steradian 
within a given redshift bin. As we determine bins with the same shot-noise, we simply have 
M = where the integral spans over the redshift range of the survey. 

To model the theoretical error , we follow the approach outlined in [62]. First, 
we take into account non-linear corrections to power spectra by rescaling all linear transfer 
functions by the fit obtained by the Haloht collaboration [52]. To do this, we have modified 
the Class code, where this rescaling for the matter power spectrum including corrections 
adapted to scenarios with massive neutrinos is already implemented [72, 73]. We extrapolate 
the Haloht rescaling, obtained only for the density term, to all the transfer functions entering 
in eq. (2.7). This gives a rough approximation for the angular-redshift power spectra C'^'^ on 
non-linear scales which is good enough for our Fisher matrix forecasts on error contours. The 
error power spectra are obtained by taking the absolute value of power spectra computed 
by multiplying all the rescaled transfer functions by the square root of 


a{k, z) 


ln[l + k/k^L{z)] , 

l + ln[l + k/k^Uz)r'"' 


(D.3) 


where A:nl('2 ) is the redshift-dependent non-linear wavenumber determined by the Haloht 
algorithm, and fth gives the error percentage on non-linear scales. Even though the claimed 
Haloht accuracy is ~ 5%, we set a more conservative fth = 10% given that the C^’s are 
computed by extrapolating the non-linear rescaling to all transfer functions (not only for the 
density). Thus, we assume that the transfer functions are affected by a fth error on non-linear 
scales k > khfhiz)- 

The Fisher matrix is then given by 


Fa0 = Yl 


dCf acf 
dXa dXp 





(D.4) 


where the Aq indicate cosmological parameters. We sum all multipoles up to £max = 2000, as 
errors on non-linear scales are taken into account via . The second sum is over the matrix 
indices (ij) with i < j and (pq) with p < q which run from 1 to Nhin when considering all bin 
auto- and cross-correlations. The hat on the covariance matrix indicates that it should be 
hrst reduced to the correlations (ij), {pq) which are considered in the sum and then inverted. 
E.g., when neglecting the signal from bin cross-correlations only i = j and p = q are kept in 
the sum, and the Cov matrix of dimension [A'bin(Abin + l)/2] x [Abin(Abin + l)/2] (for fixed 
i) is hrst reduced to a matrix Cov of dimension Abin x Abin and then inverted. 
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When constraining a subset of cosmological parameters we marginalize over the remain¬ 
ing ones. Hence, confidence levels are determined by [5] 


Ax" = E 

a,/3 


Xa) 



(A/3-A/3) , 


(D.5) 


where the hat again indicates the sub-matrix of F~^ obtained by considering only lines and 
columns corresponding to the sub-set of parameters under considerations. Fiducial values 
are indicated by an overbar Aq. The sum runs over the parameters being constrained. For 

1- dimensional and 2-dimensional contours, l-cr (68.3%) confidence levels are determined by 

Ax^ = 1 and 2.30, respectively. Note that the marginalized 1-dimensional error on the 
parameter Aq, is simply given by cTq, = while we compute the 2-dimensional 

contours shown in Figs. 6 and 8 by solving eq. (D.5). 

Finally, numerical derivatives in eq. (D.4) are computed via the five-point stencil algo¬ 
rithm, see Table 25.2 of [74]: 

... + 2K) + SCyiXa + /la) - 80^ (X^ - ha) + (Xa - 2ha) 

dXa ^ I2h ■ K ■ ) 

Here we indicate the dependence of the spectra Cf on cosmological parameters Aq. Besides 
being affected by small numerical errors scaling with steps and fifth derivatives as ^^^(c)) 

where c G [Aq — 2ha, Aq -|- 2 /iq], the algorithm is less sensitive to the ha step choice than the 

2- point derivatives. We choose /iq of the same order as the 1-dimensional errors (in this case 
obtained without marginalization) and we have verified that these are stable when further 
reducing the steps. 
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